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INTRODUCTION 


It has been documented that AE signals propagate in thin plates as 
extensional and flexural plate modes. This was demonstrated using simu- 
lated AE sources (pencil lead breaks) by Gorman [1] on thin aluminum and 
gr/ep composite plates and by Gorman and Prosser [2] on thin aluminum 
plates. A typical signal from a pencil lead break source which identifies 
these two modes is shown in Figure 1 . AE signals from transverse matrix 
cracking sources in gr/ep composite plates were also shown to propagate 
as plate modes by Gorman and Ziola [3] . Smith [4] showed that crack growth 
events in thin aluminum plates under spectrum fatigue loading produced 
signals that propagated as plate modes. Additionally, Prosser et al . [5] 
showed that AE signals propagated as plate modes in a thin walled composite 
tube . 

This fact has important implications for the interpretation of AE data 
in thin plates, shells, and tubes. First, it has been demonstrated by Gor- 
man and Prosser [2] that the source orientation can be determined by anal- 
ysis of the plate mode amplitudes. Smith [4] pointed out how this could 
be used to discriminate real AE events from extraneous noise events. Sec- 
ond, Gorman [1] discussed how erroneous source location could be obtained 
using conventional first threshold crossing or peak arrival techniques 
because of the presence of plate modes which propagate with different and 
dispersive velocities. Such source location errors were substantiated by 
Ziola and Gorman [6] and an alternative method for source location based 
on cross-correlation of the flexural mode waves was demonstrated. 

In order to apply the flexural mode cross-correlation source location 
technique, accurate knowledge of the dispersion of the flexural plate mode 
is needed. Furthermore, it has been pointed out by Mai et al . [7], Veidt 
and Sayir [8], and Dean [9] that measurements of flexural mode dispersion 
might be useful in determining the elastic constants of composite plates. 
In this research, measurements were made of the flexural mode dispersion 
on four different composite laminates. The ply layups for these laminates 
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Figure 1 Typical AE signal in thin plate identifying extension- 
al and flexural plate modes. 

were [0 16 ] , [0 4 , 90 4 ] s , [0, 90] 4s , and [0, 45, -45, 90] 2s • A Fourier phase 

technique was used to measure the dispersion up to a frequency of 160 kHz. 
on signals generated by pencil lead breaks (Hsu-Neilsen sources) . 

For future applications of this source location technique, it would 
be more convenient to analytically predict the dispersion behavior for a 
given laminate with known material properties rather than make experimen- 
tal measurements. The ability of existing plate theories for predicting 
flexural mode dispersion was investigated. The experimentally measured 
dispersion curves were compared with theoretical predictions based on 
classical plate theory (CPT) using elastic moduli calculated from lami- 
nated plate theory. The lack of agreement between theory and experiment 
at the higher frequencies demonstrated the limitations of CPT for compos- 
ite materials. These are caused by neglecting the effects of shear defor- 
mation and rotatory inertia. A higher order plate theory (HOPT) which 
includes these effects was also used to predict the dispersion behavior 
of this mode. The predictions based on the HOPT were in much better agree- 
ment with the experimental measurements. 

THEORY 

In this research, two theoretical approaches were used for predicting 
the dispersion of flexural mode waves in gr/ep composite laminates. The 
first theoretical predictions were based on classical plate theory (CPT) . 
This is a widely used approximate theory for describing motion in thin 
plates where the wavelength (X) is much larger than the plate thickness 
(h) . For flexural waves, the plate is assumed to be under a state of pure 




bending in which plane sections of the plate remain plane and perpendic- 
ular to the midplane of the plate. Thus, shear deformation is not included 
in this theory. A state of plane stress is also assumed and the effects 
of rotatory inertia are neglected. A number of authors have presented CPT 
in detail including Graff [10] who derives the equation of motion for iso- 
tropic materials and Whitney [11] who includes the effects of anisotropy. 
The CPT equation of motion for an orthotropic composite laminate in the 
absence of body forces is 
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where the D i j , s are the anisotropic bending stiffness coefficients 
obtained from laminated plate theory as described by Whitney [11] or Tsai 
and Hahn [12] . In the previous equation, w is the displacement along the 
z axis or normal to the plane of the plate, x and y are orthogonal axes 
in the plane of the plate, and p is the density. 

The dispersion behavior for the flexural mode using CPT is obtained 
by substituting the displacement for a plane wave propagating in an arbi- 
trary direction into the equation of motion. This displacement is of the 
form 


w = A 0 e i(cot - klxX - kl >- y) Eg. 2 

where A 0 is the amplitude, oo is the angular frequency, l x and l y are the 
direction cosines of the direction of propagation, and k is the wavenum- 
ber. After substitution and reduction of terms, the resulting CPT dis- 
persion relation is found to be 
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where c f is the velocity of the flexural mode and 
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Thus, in CPT, the velocity is dependent on the direction of propagation 
and increases as the square root of the frequency without limit. 

The second theory used to predict the dispersion was a higher order 
plate theory (HOPT) which includes the effects of shear deformation and 
rotatory inertia. This theory was put forth by Tang et al . [13] following 
earlier work by Yang et al [14] which was an extension of work by Mindlin 
[15]. The dispersion behavior for a symmetric orthotropic laminate pre- 
dicted by this theory is obtained when the determinant of the following 
matrix of coefficients is set equal to zero 
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In Eq. 16, the Kq are shear correction factors which were determined to 
yield the best aqreement with three dimensional elasticity theory when 
K i 2 =5/6. The subscript 1 refers to the l'th layer of the laminate and the 



Qj_j are the stiffnesses for the l'th layer. Solving the determinant for 
the wavenumber as a function of co yields a cubic in k 2 . Only the root 
which approaches zero as the frequency approaches zero is the correct 
root. Once k as a function of co is known, the phase velocity is determined 
as a function of frequency using the relation 

= 7 - . Eq. 17 
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EXPERIMENT 


The composite laminates used in this study were made of AS4/3502 
graphite/epoxy. All four laminates consisted of sixteen plys and had a 
nominal thickness of 2.26 mm. The dimensions were 0.508 m. along the x 
direction (0 degree ply direction) and 0.381 m. along the y direction. 
Measurements were made along the 0 degree (x direction) , 45, and 90 degree 
directions for all four laminates. The nominal lamina properties for this 
material as obtained from the manufacturer are given in Table 1 . These 
values were used in the laminated plate theory calculations to obtain the 
bending stiffness coefficients needed for the CPT and HOPT dispersion cal- 
culations . 


Lamina thickness = 1.413 X 10 4 m. 

_ o 

Density = 1550 kg/m Fiber volume = 60% 

Q xx = 145.5 GPa . Q xy = 2.91 GPa Q yy = 9.69 GPa Q ss = 5.97 GPa 


Table 1 Lamina properties of AS4/3502 graphite epoxy. 

A Fourier phase technique was used for the measurement of the flexural 
mode dispersion. This technique has been described by a number of authors 
including Sachse and Pao [16], Pao and Sachse [17], Veidt and Sayir [ 8 ], 
Dean [9], and Alleyne and Cawley [18]. In this technique, the elastic 
wave is detected at two different distances away from the source of the 
wave along the direction of propagation of interest. The phase (cp) of 
the wave at each position at a given frequency (f) is determined by per- 
forming a Fourier Transform on the signals and computing the phase. The 
phase must be unwrapped to remove the 2njt uncertainty. The phase differ- 
ence (Acp) over the distance between the two transducers (Ax) is then com- 
puted for each frequency. The wave number and velocity are then calculated 
at each frequency by 
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The experimental setup used for these experiments is shown in Figure 
2. The two receiving sensors were Panametrics 3.5 MHz broad band ultra- 



Figure 2 Experimental setup for flexural velocity measurements 
in composite plates. 


sonic transducers. These sensors have been shown by Prosser [19] to pro- 
vide flat frequency, displacement sensitivity response to these low 
frequency plate waves. The trigger sensor was a Physical Acoustics Cor- 
poration (PAC) model R15. The preamplifiers (PAC model 1220A) were set 
at 40 dB amplification with no filtering. The source was a pencil lead 
break repeated with the transducers at separations of 1.91, 2.54, 3.18, 
3.81, and 4.45 cm. An average velocity and standard deviation for the 
five different measurements was computed. The source and receivers were 
kept as nearly in the center of the plate as possible to minimize reflec- 
tions. The waveforms were digitized at a sampling frequency of 1 MHz with 
a LeCroy 6810 transient recorder and then transferred to a personal com- 
puter for processing. 

Prior to computing the FFT to determine the phase, the higher fre- 
quency extensional mode and the reflections arriving later in the flexural 
mode were zeroed out in the computer. Previous Fourier analysis of the 
flexural mode signals when digitized at much higher sampling frequencies 
(100 MHz) showed that the maximum frequency component in the flexural mode 
was about 200 kHz. Thus, aliasing was not a concern even at the low sam- 
pling frequency of 1 MHz. 


RESULTS AND DISCUSSION 


The average measured velocities for the 0, 45, and 90 degree direc- 
tions in the [0 16 ] graphite epoxy plate are plotted in Figure 3 to Figure 
5 with the standard deviation of the measured values indicated by error 
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Figure 3 Measured and theoretical flexural dispersion for 0 de- 
gree propagation in [0 16 ] graphite/epoxy plate. 

bars. The predicted velocity dispersion curves for CPT and HOPT are also 
shown in these plots. The agreement between measurement and HOPT is excel- 
lent for the 90 degree propagation direction. For propagation at 45 and 0 
degrees, the measured values are consistently less than those predicted 
by HOPT. The discrepancy in the HOPT predicted and measured velocities 
at 0 and 45 degrees was believed to be caused by differences in actual 
material properties from the nominal ones used in the calculations. This 
is common for these materials and is due to fiber volume variations, cure 
processing variations, and variations in resin chemistry. These discrep- 
ancies were consistent with differences in predicted and measured exten- 
sional velocities for these same plates reported by Prosser [19]. 

The effect of shear and rotatory inertia is clear when the CPT and the 
HOPT are compared in these plots. CPT and HOPT are in agreement at very 
low frequencies in all cases where the approximations of CPT are valid. 
The discrepancy between the two increases with increasing frequency as the 
velocity predicted by CPT increases without bound. 

It is also apparent that the difference between HOPT and CPT is much 
greater for the 0 and 45 degree directions than for the 90 degree direc- 
tion. This is expected since the shear modulus is much smaller in com- 
parison to the Young's modulus in those directions. CPT is based on the 
assumption of no shear deformation which implies an infinite shear modu- 
lus. Thus, better agreement is provided by CPT when the ratio of the shear 
modulus to the Young's modulus is larger. 





Velocity (m/s) Velocity (m/s) 






A plot of the measured velocities and the HOPT predictions for the 
[0, 90] 4 g plate is shown in Figure 6. In order to better view this com- 
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Figure 6 Measured and HOPT theoretical flexural dispersion for 
0, 45, and 90 degree propagation directions in 
[ 0 , 90] 4 s graphite/epoxy plate. 


plicated graph, the CPT predictions and the experimental uncertainties are 
not shown in the plot. In this plate, the measured velocities were less 
than theoretical predictions for all three directions of propagation. 
However, the measured and theoretical velocities occurred in the same 
order with the 0 degree velocity being the largest, followed next by the 
90 degree velocity, and with the 45 degree velocity the smallest. The 
results for the other two laminates were similar with the measured veloc- 
ities consistently less than predicted by HOPT. This again seems to indi- 
cate that the actual material properties are less than the nominal 
properties used in the theoretical calculations. 

In summary, a Fourier phase technique was used to measure the disper- 
sion of the flexural plate mode in four composite laminates. CPT was shown 
to be of limited value in predicting the dispersion of this mode because 
it assumes that the effects of shear deformation and rotatory inertia are 
negligible. HOPT, which includes these effects, gave much better agree- 
ment with the measured values. However, there was a consistent discrep- 
ancy between theory and experiment believed to be due to variations in 
actual material properties from those used in the calculations. Thus, 
HOPT looks promising for predicting the dispersion behavior of the flex- 
ural mode for use in the cross-correlation source location technique. 




However, accurate material properties are needed. If these are not known, 
an experimental technique for measuring this dispersion has been pre- 
sented. 
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